Water waves over a strongly undulating bottom 
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Two-dimensional free-surface potential flows of an ideal fluid over a strongly inhomogeneous bot- 
tom are investigated with the help of conformal mappings. Weakly-nonlinear and exact nonlinear 
equations of motion are derived by the variational method for arbitrary seabed shape parameter- 
ized by an analytical function. As applications of this theory, band structure of linear waves over 
periodic bottoms is calculated and evolution of a strong solitary wave running from a deep region 
to a shallow region is numerically simulated. 
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I. INTRODUCTION 

The classical problem of water waves over a variable 
seabed has attracted much attention (se e El B IM II, IE 

lillia 13 m Ha |21, 22., ^, ^ and references therein) . 
There are some significant differences in this interesting 
and practically important problem, as compared to the 
theory of waves on a deep water or in canals with a flat 
horizontal bottom. In situations where the fluid depth 
is less or of the same order as a typical length of sur- 
face wave, inhomogeneity of the bottom is a reason for 
linear and nonlinear wave scattering and transformation, 
and it strongly affects wave propagation. These phenom- 
ena occur so widely that one can meet them almost ev- 
erywhere, although with different scales. Examples of 
strongly nonlinear dynamics are ocean waves running on 
a beach, or motion of disturbed water in a puddle after 
a car. Among linear effects due to bottom topography 
is existence of special edge-localized waves discovered by 
Stokes |i3, m, Ea I2OI , that propagate along the shore 
line of a beach. Over an axially symmetric underwater 
hill, quasi-localized wave modes with non-zero angular 
momentum can exist, similar to long-life-time states of a 
quantum particle confined by a potential barrier of a fi- 
nite width [a IsJ • It is necessary to say that underwater 
obstacles of definite shapes and sizes can serve as waveg- 
uides (a narrow and long underwater crest) or as lenses 
(an oblong underwater hill oriented crosswise to the wave 
propagation). A qualitative explanation for all the lin- 
ear effects is simple. Indeed, let rj_ be the coordinate in 
the horizontal plane, H{r±) the depth corresponding to 
quiet surface. Then, looking at the well known dispersion 
relation for small-amplitude gravitational surface waves. 



uj{K,H) = ^/gK tanhiKH) 



(1) 



(where to is the frequency, K is the absolute value of 
the wave vector, g is the gravitational acceleration), one 
can see that the local refraction index n{uj, r±) increases 



as the depth H{r±) decreases, in accordance with the 
formulas 






gK{uj,H{v^)) 



K{cj,H = 00) 



>1, 



(2) 



dK{uj,H) 
'dH 



<0, 
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where the function K{uj,H) is determined by Eq. ^. 
Therefore, as in the conventional light optics, here an 
oblique wave changes its direction of propagation when 
meets gradient of n. Also, the total internal reflection 
is possible in propagation from smaller depth to larger 
depth. 

Besides observing such natural phenomena, a set of 
laboratory experiments has been carried out to investi- 
gate various aspects of the given problem in more ideal- 
ized and controlled conditions than are achieved in na- 
ture j9,j 12, 13; 14, 19, 20, 27]. In particular, waves 
over locally periodic bottoms were studied experimen- 
tally [^ llj, 113, m, I23 , and such a general for periodic 
media effect was observed as the Bragg resonances and 
the corresponding band structure with gaps in wave spec- 
trum. It is worth to say that in natural conditions quasi- 
periodic sand bars occur quite often. 

In general, a qualitative picture of the mentioned phe- 
nomena is clear. As concerning the quantitative side of 
the mathematical theory of waves over a variable bot- 
tom, here not everything that necessary has been done, 
because practically all developed up to now analytical 
models and methods are related to the limit cases where 
the fluid is considered as ideal, and the slope of the bot- 
tom is small (or amplitude of the bottom undulations is 
small). For the general three-dimensional (3D) Hamilto- 
nian theory of water waves, such restriction seems to be 
unavoidable even in considering the most simple, irrota- 
tional flows when the state of the system is described by a 
minimal set of functions, namely by a pair of canonically 
conjugated quantities as the deviation ?7(r^, i) of the free 
surface from the horizontal plane and the boundary value 
'(/'(r_L, t) of the velocity potential 32, 33] . A technical dif- 
ficulty exists here that, when working in 3D space, it is 
impossible to represent in convenient and compact form 



the kinetic energy functional JC{r], tp} which is part of 
the Hamiltonian of the system. Small values of the bot- 
tom slope and of the free surface slope make possible 
expansion of the Hamiltonian to asymptotic series and 
subsequent application of various variants of the pertur- 
bation theory. In such traditional approach, an inhomo- 
geneous bottom does not allow to write in exact form 
even linearized equations, not speaking about nonlinear 
corrections. 

There are more favorable conditions for progress in 
theory of 2D potential ideal flows with a free bound- 
ary, and the reason for this is the possibility to em- 
ploy such powerful mathematical tools as analytical 
functions and the corresponding conformal mappings. 
Time-dependent conformal mappings were successfully 
used for studying strongly nonlinear 2D wave dynam- 
ics on deep water and over straight horizontal bottom 
Hi m Hi El S m ES El- In the cited works the 
region occupied by resting fluid (the lower half-plane or 
a horizontal stripe) was mapped onto the region with dis- 
turbed free boundary, and the real axis was transformed 
into moving boundary. Such a conformal "straightening" 
of free surface has provided a compact representation for 
the Hamiltonian, derivation of exact equations of mo- 
tion, and possibility for precise numerical simulations of 
the system evolution. 

The purpose of this work is to study the effect of a 
strongly undulating bottom on 2D ideal potential flows 
with a free surface. Here conformal mappings are used 
as well, and this is done in two variants. In the first, 
"moderate" variant (Sec. 2), a fixed conformal map- 
ping "straightens" the bottom, but not the free bound- 
ary. More exactly: instead of the Cartesian coordinates 
X and y (with y-axis up-directed), curvilinear coordinates 
u and V are introduced, and the change of coordinates is 
performed with the help of an analytical function z{w) 
which maps the stripe — 1 < Imw < onto the region 
between the horizontal line y — and the inhomogeneous 
bottom y — —H{x). In this case x -\- iy = z{u -\- iv), the 
horizontal line y = corresponds to w = 0, and on the 
bottom V = —1. The bottom may have arbitrary large 
slope and even impending pieces where the dependence 
H{x) is multi- valued, as shown in Fig^ The shape of free 
surface will be described by a function v = V{u,t). The 
Lagrangian for weakly-nonlinear waves is represented as 
an integral series in powers of the dynamical variables 
V(u,t) and ip{u,t), with coefficients explicitly depend- 
ing on the spatial coordinate u. In the small-amplitude 
limit, the wave dynamics is governed by linear integral- 
differential equations. It is using the conformal variables 
u and V, that allows us to obtain these equations in 
exact form, contrary to the traditional approach where 
even linearized equations can be obtained only approx- 
imately by expansion in the small parameter, the slope 
of the bottom. The definition "moderate" for this vari- 
ant emphasizes that straightening of the bottom without 
straightening the free boundary is able to provide not 
more than a weakly-nonlinear theory. Nevertheless, such 




FIG. 1: Left: periodic shape of the bottom (lower thick line), 
levels of constant v = —0.9, —0.8, ..., -1-0.6 (thin lines), and 
(schematically) free boundary (thick line near y — 0). Right: 
the corresponding band structure of the spectrum of linear 
waves. In this example z{'w)/h = w + {2e/a) sin(Q:«;)/(l -|- 
bcos{aw)), with the parameters a — 2-k/3, e — —0.17, b — 
0.16. 



a theory seems to be helpful and applicable in many prac- 
tical cases when wave amplitude is small. The results of 
this part of the work are the derivation of the Hamil- 
tonian functional for weakly-nonlinear potential surface 
waves in canals having arbitrary bottom shape, as well as 
calculations for band structure of spectrum for a number 
of periodic bottom profiles. As an example how to treat 
the linearized equations, also the problem is considered of 
wave reflection on a smooth "step" - changing the depth 
from hi to /i2. 

The other variant of using the conformal mappings 
may be called "radical" in the sense it is valid for ar- 
bitrary shape of the bottom and for arbitrary shape of 
the free surface. It is an exact combined theory where a 
time-dependent conformal mapping straightens both the 
bottom and the free boundary (Sec. 3). Such a mapping 
can be represented as the result of two mappings: x+iy — 
z{C{'w, t)), where the flrst function C(w', t) maps the hor- 
izontal stripe — 1 < w < onto the region D(^{t) with 
the straight lower boundary (Im(^ — —1) and with a per- 
turbed upper boundary, after that the time-independent 
function z(Q maps the half-plane Im(^ > — 1 onto the re- 
gion y > —H{x) in the physical plane bounded from be- 



low by the bottom. The shape of the free surface will be 
described by the formula X + iY — Z{u,t) = z{({u,t)). 
However, it appears that exact nonlinear equations for 
Z(u, t) in the inhomogeneous case have the same form as 
the known equations for waves over a horizontal bottom 
|38| . but with different analyticity requirements imposed 
on the solutions. Numerical solutions obtained by the 
spectral method are presented that describe a running 
and breaking wave (Sec. 4). 



II. WEAKLY-NONLINEAR THEORY 

So, suppose we know the analytical function z{w) = 
x(u, v) + iy{u, v) which maps the horizontal stripe —1 < 
Imw < onto the region occupied by the fluid at rest, 
and this function takes real values on the real axis: 
z{u) = x{u, 0). The velocity field is irrotational, and the 
velocity potential Lp{u, v) satisfies the Laplace equation 
fuu + ^vv = in the flow region —l<v< V{u, t), with 
the boundary conditions ipv\y^-i = 0, (p\v=v{u) — i^i"^)- 
Due to conformal invariance of the Laplace equation in 
2D-space, hence equation Lp^x + '^yy = is satisfied as 
well, with no-penetration boundary condition on the bot- 
tom: dip/dn\y^_H(x) = 0. Let us now take into account 
the fact that the Lagrangian functional for potential sur- 
face waves has the following structure |32. l33l i38.] : 



C 



tprjclx — Ti. 



iPiYtXu - YuXt)du - n, (3) 



where Y{u,t) — y{u,V{u,t)), X{u,t) = x{u,V{u,t)), 
and the Hamiltonian functional Ti is the total energy of 
the system - sum of the kinetic energy and the poten- 
tial energy in gravitational field (in this paper we neglect 
surface tension effects, though they can be easily incorpo- 
rated by adding to the Hamiltonian the surface energy). 
In our variables 



1 



V{u) 



n = -J du J {ipi + ipi)dv 
-1 

^ f y\u,V{u))-^x{u,V{u))du. (4) 



where x'{u) = z'(u-\-Oi) = x^iu, 0), and the equality yy — 
Xu has been taken into account in the first integral in 
r.h.s. The expansion for the kinetic energy (calculation of 
the functionals /C*^^^ and /C^'^^) is performed in a standard 
manner 33. 38J and gives 



^(2){^} 



^[k tank k]^du, 



(6) 



/C(3){^,T/} 



[4>l-{[ktaTihk]iPf]Vdu. (7) 



Here the linear Hermitian operator [k tanh k] has been 
introduced, acting as 



-t-oc 



[fc tanh fc] ■!/;(«) 



P.V. 



ipu{u)du 



2sinh[(7r/2)('u--'u)]' 



(8) 



In Fourier-representation this operator simply multi- 
plies the Fourier-harmonics tpk — J ^{u) exp{—iku)du by 
A:tanhfc. Quadratic on V part of the potential energy is 



V^^\V}=^-jv^x'\u)du. 



(9) 



It is convenient to deal with the function £,{u,t) = 
V{u,t)x''^{u) canonically conjugated to ilj{u,t), and 
write the corresponding up-to-third-order Hamiltonian in 
terms of ^ and ip: 



7Y{C,^} 



i/#tanhfc]^..+ f 



e 



x'(u) 



du 



[ipl~ {[kta.-Dhk]%l)f]S, 



du. 



(10) 



Physically, this asymptotic expansion of the Hamiltonian 
is on a small parameter - the slope of the free surface (see 
J33| for more comments and references). The weakly- 
nonlinear equations of motion have the standard Hamil- 
tonian structure 

Ct = ^ = [kanhfclV - — (^^ 
dip du \x'^{u) 



[k tanh k] 



^[kta,nhk]tp 
x'^u) 



(11) 



This system has the obvious stable equilibrium ip — 0, 
V — 0, hence one may consider weak oscillations near this 
equilibrium state. In a standard way (see, for instance 
|33]), let us expand the Lagrangian © in powers of the 
dynamical variables ip and V. It is clear that due to 
the symmetry principle the expansion for j/(m, v) contains 
only the odd powers of u, while the expansion for a;(M, v) 
contains only the even powers of v. Therefore up to the 
third order in powers of ip and V the Lagrangian ^ is 
equal to 



c = I ipvtx''^{u)du - a:(^){v} - /c^^^iv-, V} - V^'^'^iV}, 

(5) 



V't^ 



^ 



sn 



5£, "x'iu) 



[V^g-([fctanhfc]Vj)^ 
2x'^{u) 



(12) 



If \x" [u) / x' {u)\ <C 1, then x'{u) is approximately equal 
to the equilibrium depth H{u). For long waves over a 
such slowly varying bottom, only Fourier-harmonics ipk 
and S,k with small k are excited, so in this case the Hamil- 
tonian (|10|l can be simplified to the local form 



Wi = 



2 



^: 



2iP. 



2 
uuu 



ge 



ai^l - 1 



^ 1 

uu\ 



6 



I I IILLLL 1^ -v ^ _i 3 L ' -(i. ' act J 7 

^ ~l5~ 2^^ ^ 2x'Hu) \ ' 

(13) 

which is suitable for consideration of such phenomena as 
interaction of solitons with the bottom topography. 



A. Linearized equations 

Now let us consider the linearized system 



ft — [fc tanhfcjV', 



-V't =9 — 



x'[u) 



(14) 



For a monochromatic wave (f , V' oc exp(— icji)) Eqs. H14|) 
are reduced to the single integral equation 



— a;'(w)-fctanhfc ■0c^(m) = O. (15) 

9 



In the low-frequency limit this equation can be consider- 
ably simplified. A variant of simplification is to introduce 
a new function / by the equality A; tanhfc?/;w(u) = —fuu- 
Then we obtain the equation 



— a;'(u)fc coth k + (d/duf ) / = 0. (16) 

9 



The low-frequency limit corresponds to long wave- 
lengths, when fccothfc w l + k'^/3 = 1~ {l/3){d/du)'^, so 
we have to deal with the second-order differential equa- 
tion 



fuu{u) + 



-x'{u) 



^x' (u) 



f{u) = 0. 



(17) 



where uj'^x'{u)/g should be small (only in this case the 
wave length is indeed effectively long; remember that 
x'{u) is of the same order as the depth). 

Higher-order approximations to equation H15|) can be 
derived in a similar manner, for instance by change 
i^uj{u) = [coshfc]/(u) and subsequent expanding [coshfc] 
and [/csinhfc] in powers of k^ — ~{d/duY. 

As an explicit example of using Ea. p7|l . we consider 
reflection of a long wave from a step-shaped bottom in- 
homogeneity described by the function 

z{w) ^ hiw + ^^' ~ ^'^ ln(l + e""'), (18) 

a 

where /ii > ft.2 > 0, < a ^ tt. If frequency of the 
wave is small, uj ^ \/g/hi, then equation (|17|l may be 
applied. Calculating the derivative 



1 



Z'{w) =/ll + (/l2-/ll)- 

we have for f{u) the equation 



hie 



(19) 






where 



/ii- 



hi 



3g 



h2 = 



1 - 



C = 



LP hi 

3g 



3ff 



1- 



3s 



A general solution for equation H20() is known [31|. In 
particular, the reflection coefficient is given by the ex- 
pression 



RiLu) = 



sinh[^{Vhi-Vh, 



,sinh[i^(V/ii + 




(22) 



B. 



Periodic bottom: The band structure of the 
spectrum 



Interesting phenomena occur if shape of the bottom is 
periodic: 

z'{w) ~ hy ^ a„ exp(mo;w), a-n = dn- (23) 

n 

Here /i is a dimensional parameter, a„ are some complex 
Fourier-coefficients. Obviously, x'{u) = z'{u) > and 
|a„| decay rapidly at large |n|, since z'(w) does not have 
any singularities at — 1 < Imw < 1. The equation H15|l 
for eigen-functions '4'\{u) (where A = uP^hj g) now has the 
form 

A y ani2xp{inau) tp{u) — [kta,nh k]'ip{u) = 0, (24) 



or in Fourier-representation 

A^a„Vfc-nQ = fctanhA: V'fe- 

n 

For convenience let us denote 

F^ = av tanh(ai^), $,y = ^'av 

Now we have the infinite chain of linear equations 



(25) 



(26) 



(27) 



(21) 



where "^ ^^ and ^^3 interact if the difference between vi 
and V2 is an integer number. Let us fix some v. Non- 
trivial solutions of the system H27|l exist only at defi- 
nite values A = Xm{v), where m = 1,2,3,.... It is 
necessary to note that the functions \m{i^) are periodic: 
\m{v -I- 1) = \m{v), and even: \m{-v) = >^m{v)- This 
determines the band structure of the spectrum with fre- 
quency gaps (see Figs lll2|) . For numerical computing 
Am [v) it is necessary to cut the infinite chain H27(l at some 
large but finite length, thus considering only v between 
—N and N . Practically N should be several times larger 
than the index m of Am . Numerical results for vAm 
shown in the figures [TO were obtained with the help of 
the mathematical package Maple 8 taking A^ = 10. 

FigHlshows that in some cases even for strongly undu- 
lating bottom the coefficients a„ with n > 1 can be still 
small (fli = e = 0.2 <C 1). In these cases it is easy to cal- 
culate analytically in the main approximation positions 




FIG. 2: Shape of the bottom, levels of constant v, and the 
band structure for z{w)/h = w + (2e/Qf) sin(aw), with a — 
7r/3, e = 0.2. 



of the gaps. For example, let us consider the bottom pro- 
file as in Fig|21 where x'{u) — h{l+2€ cos{au)) . The gaps 
in spectrum correspond to integer or half-integer i^'s. It 
is important that at these values of v, solutions of the 
linear chain H27|l possess definite parity, in the sense that 
'i'-v = ±^iy. This allows us in gap calculation consider 
only positive v. Let us first consider half-integer z/'s and 
the corresponding half-infinite chain 



(A - i^l/2)*l/2 + Ae(±*i/2 + *3/2) 

(A - ^3/2)^3/2 + Ae(*i/2 

(A - F5/2)*5/2 + Ae(*3/2 



*5/2) 
*7/2) 



0, (28) 
0, (29) 
0, (30) 



Obviously, the even and odd cases result in different A's, 
and it is this difference that determines the gaps in spec- 
trum. For main-order calculation of the first and third 



gaps, we cut this chain: \1/ 



7/2 



0, ^ 



9/2 



— 0, and so 



on. Now we have to solve the equation for zeros of the 
determinant 3x3 



[{A(l ± e) - F1/2XA - F3/2) - A^e1(A - F5/2) 
-A2e2[A(l±e)-Fi/2]=0. 



(31) 



First we take A = i^i/2-t-Ai, where Ai is a small quantity 
of the order e. In the main order Ai ± e_F\ /2 ~ 0, and 
this gives us the first gap: i^i/2(l — e) < A < -^1/2(1 + e)- 



For the third gap we write A = i^3/2 + A3, where A3 
is of order e^ . The equation for A3 with the third order 
accuracy is 



[(F3/2(l ±e)- i^i/2)A3 - 6^^3^/21(^3/2 - F, 



'5/2) 



-£'^3/2(^3/2(1 ±e)-Fi/2) = 0. (32) 



From here we find 



A3 = e'F. 



3/2 



(^3/2-^^5/2) (F3/2(l±e)-Fi/2) 



where we may keep only the second- and third-order 
terms. This gives us the position of the third gap 
AL'UA<Af: 



Al^) = K 



3/2 



±- 



*= -^3/2 



3/2 



(-P3/2 - ^5/2) 



■(i^3/2~i^l/2)' 



(^3/2 - i^l/2) 

(33) 



Analogously, the gaps at integer v's can be considered. 
These are determined by the system 

(A-Fo)*o + Ae(±*i + *i) = 0, (34) 

(A-Fi)*i + Ae(*o + *2) = 0, (35) 

(A-F2)*2 + Ae(*i + *3) = 0, (36) 



For instance, position of the second gap in second order 
is given by the formulas 



aL') = f^- 



F2 — Fi' 



\f =Fl{\ + 2e^)- 



e'Fl 



Fo 



Fi 
(37) 



III. EXACT THEORY 



In exact nonlinear theory, the shape of free bound- 
ary is given in parametric form by a compound func- 
tion z(C(u,t)), where z(C) is a known function com- 
pletely determined by the bottom shape [for example, 

z{() — /i(\/(C + *)^ ~ {b/h)'^ — i) corresponds to a narrow 
vertical barrier of the height & at a; = on the straight 
horizontal bottom with the depth y = —h]. The un- 
known function C(z/;, i) should be analytical in the stripe 
— 1 < Im w < and the combination [(^{u—i,t)+i] should 
take real values. These conditions relate the real and the 
imaginary parts of C{u,t) at the real axis |3q |: 

C{u,t)=u+{l + iR)p{u,t), (38) 

where p(u, t) is a real function, and the linear 
anti-Hermitian operator R is i tanh k in Fourier- 
representation. In u-representation 



-I-CX3 



Rp{u,t) =P.V. 



p{u,t)du 



2sinh[(7r/2)(w-u)]' 



(39) 



The inverse operator R ^ = T 

-\-oo 

fp{u,t) = P.V. 



: —i coth k acts as 
p{u, t)du 



1 — exp[7r(u — u)] 



(40) 



Note that the previously considered operator [/ctanhfc] 
is —Rdu- The kinetic energy functional is now exactly 
equal to the expression at the r. h. s. of Eq.®. The 
Lagrangian for ^{u, t) and C,{u, t) is given by the formula 



C, 



F (01 -. ) Hu 



A 



^pRip^du 



2i 



{Imz(C)}2Re(z'(C)C«)rf" 



2i 



-R 



c + c 



du, 



(41) 



where the (real) Lagrangian indefinite multiplier A(u,i) 
has been introduced in order to take into account the 
analytical properties of the function C, given by Ea. (|38|l . 
From the above Lagrangian one can obtain the equations 
of motion. Variation of the action J Ccxactdt by (5^(u, t) 
gives the equation 



\z'{C)\'{Ut-Ut)/{2i)^~Ri^u, 



(42) 



which can be easily transformed to the form (compare 
with HI) 



Ct^-UT + ^) 



Ripu 



WiCWKul 



(43) 



The variation of the action by dC,{u,t) results after sim- 
plifying in the equation 

k'(C)l' {V-tCn - VnCt +3CJmz(C)} - (1 + lR)A - 0. 

(44) 
Since the product Cm(1 + i-R)A has the same analytical 
properties as both C„ and (1 + iR)A, we can multiply 
Ea. (|H|l by Cm and write 

|2'(C)|' {[A + <?Imz(C)]|CnP - V'nOC.} - (1 + *^)A = 0, 

(45) 
where A is another real function. The imaginary part of 
the above equation together with Ea. (|42ll result in 



A=-f[V'ni?VJ- 



(46) 



Using this equality, we can reduce the real part of Eq. ||1S|) 
to the form 



ipt +glT[iz{C) = -TpuT 



R^u 



WiCWKu 



flTpuR-ipu] 



(47) 
which is the Bernoulli equation in the conformal vari- 
ables. Exact equations (|43|l and H47|) [with given analyti- 
cal function z{Q and with the condition l|38|) ] completely 



determine the evolution of gravitational surface waves 
over the undulating bottom parameterized by a real pa- 
rameter r as Xi,{r) +iYi,{r) — z(C)|(^=r-i- 

Equations 143|) and (|47|l can be represented in another 

form by using the identity 2T[^puRipu] = ipu — {Ripu)'^ 
|38| and introducing the complex potential 



$(u,i) = {I + iR)^p{u,t) 



(48) 



(which is analytically continued to the stripe —1 < 
Imw < 0): 



Ct = -UT + i) 

$t - -$„(r-hz) 

~{l + iR) 



Im$„ 



^'(OPICuP 

Ini $„ 



k'(C)PIC«P 
2k'(C)PIC«P 



5lniz(C) 



(49) 



(50) 



A very interesting point is that one can re- write equations 
(|43|l and (|47|l without the intermediate function C,{u,t), 
but directly for z{u,t). Indeed, after multiplying Eg. jl^ 
by z'(Q we obtain the equations 



zt = -Zu{T + i) 



\Zu? 



V't -l-glmz = -i^uT 



\Zu\' 



T[lpuRll^u] 



(51) 



(52) 



that is exactly the same system as was derived in |2a| 
for a straight horizontal bottom. However, in our case 
analytical properties of the function z{w,t) are different: 



Imz(u) y^ R[Re {z{u) — u)]. 



(53) 



The only requirements for z{w,t) now are that it should 
be analytical in the stripe — 1 < Imw < and the cor- 
responding mapping should have a physical sense (no 
self- intersections are allowed). The question may arise: 
Where is the bottom shape in Eqs. (|FT|l - lf5^ ? The an- 
swer is simple: The shape of the bottom is an integral of 
motion for this system. Roughly speaking, each particu- 
lar solution of Eas. H51|l - l|52|) corresponds to a flow over a 
definite topography determined by the initial condition 
z{r — i, 0). 

Analogously, Eas. lfl^ - lfSUIl can be represented as 



Zt 



$* 



-$„(T + i) 



-- -zu{i 


+ *) 


'Im$„ 

[ \Zu? 


(54) 


"Im$„' 


~{\ + iR) 


\'^n? , , 1 




[ \ZuV \ 


2|z„P ' ^^""J 
(5 
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IV. NUMERICAL EXPERIMENTS 

A. Different forms of equations 

For numerical simulations, still other equivalent forms 
of exact equations may be useful, since numerical sta- 
bility depends dramatically on the choice of dynamical 
variables j^E 1^ Ell ^M ■ Two alternative sets of equa- 
tions were used in computations presented below. First, 
as it was pointed in [33 for the case of deep water, a good 
practical choice for the dynamical variables is A = l/z„ 
and B = $„/z„. It is easy to derive the equations of mo- 
tion for A{u,t) and B{u,t) from Eas. (|5^ - lf55|) . and they 
are very elegant (compare with 



7i-period. So, at arbitrary moment of time we will have 
the equality C(w,i) +i = (-L/27r)C* {2Tra{t)u/L,t), where 



+00 



^ ^ ^ ^' ^-^ l + exp(2ma(i)) 

771— — 00 



(61) 



with an unknown real function a(i). The unknown com- 
plex Fourier coefficients yOm(i) correspond to a real (27r- 
periodic on the variable ■&) function p{'&,t): 

-t-00 

pid.t)^ ^ p„(t)exp(zTOZ?), P-m{t) = Pm{t). 



At = -Auif + i)lm (BA) + A{f + i)9,Jm (BA), (56) 



Bt = -BuiT+t)liJi{BA)-A{l+iR) 



-^ + ^''''a 

(57) 
The variables A and B do allow stable numerical simula- 
tions for waves over varying seabed. However, analytical 
properties of A and B are not restricted by conditions 
similar to Ea. (|38|l . and therefore the shape of the bot- 
tom is preserved in this case only approximately. 

The second set of variables, that were used in numeri- 
cal experiment, consists of two complex functions: C(u, t) 
and 



/3(u,i) = ^uiu,t)/Cu{u,t), 



(58) 



both having effectively controlled analytical properties. 
With this choice, the bottom shape is preserved exactly, 
but the corresponding equations of motion are slightly 
less compact: 



(3t 



-Cn(T + i)Im 
-/3„(f + i)Im 



p 



p 



2k'(C)P 



+ gluYz{C,) 



(59) 



(60) 



It is necessary to explain here some important details 
about space-periodic solutions of the system (|55|) - (|5n|l . 
since spectral numerical methods deal with periodic func- 
tions. Such solutions exist if the function z'{C) is periodic 
with a fixed real period L, so that z{C, + L) = L + z{(). 
However, this does not imply that the functions C„(u,t) 
and P{u, t) have a fixed w-period. It would be so, but the 
linear operator T is singular at small fc, and its action 
on a constant function is not periodic in u: TC = Cu. 
Thus in the right-hand-sides of the Eas. (|59|l - H60l) we have 
non-periodic terms. Therefore C,t(u,t) and Pt{u,t) can- 
not retain a constant w-period. However, (^tiuji) and 
Pt{u,t) can be space-periodic with a time-dependent 



The linear operator R^ is diagonal in the discrete Fourier 
representation: Rq(to) — itanh(Q;7TT.). 

Analogously, f}{u,t) can be represented as (3 — 
(<7L/(27r))i/2/34^,i), where 



+ 00 



P*i^,t)^ J2 Y 



+ exp{2ma{t)) 



(l + iR„)x(^,i). 



(62) 
Equations of motion for the real functions a(t), p{'d,t)^ 
and x(^,i) follow from Eas. (|S5|l - l|Hn|l : 



a{t) 



1 



Im 



~IW) 



2vrjo 
p(^,t) = Refc:(Ta+*)Im 






M, 



x(^,i) = Re(/3:(t„ + 7;)Im 



--:7(l + iRa)a^ 



i<(c*)pc:(^) 
-p*m 

K{C*)\H'M 



2|4(C*)l^ 



+ Imz*(C,) 



(63) 
(64) 

(65) 



where z*(C*) — {2-K/L)z{—i -f L(^./{2tt)), (^;|^ — 9^C* £^nd 
so on. The linear operator T^ is regular. In the discrete 
Fourier representation it is defined as follows: 



Ta(m) = — icoth(aTO), 
= 0, 



TO^O; 

m — 0. (66) 



Two numerical experiments are briefiy reported below, 
first of them employing Eas. H56l) - (|57|l . and the second one 
employing Eqs.|(5ll)-l|nnil- Both the systems (|5S1)-(|57)) and 
(|59|I - (|6U|) are equally convenient for numerical solution by 
spectral methods, inasmuch as the multiplications can be 
performed in u-representation while the linear operators 
R and T (also the ^-differentiation) are simple in the 
Fourier representation. Efficient subroutine libraries for 
the fast Fourier transform are now available. The inte- 
gration schemes in both cases were based on the Runge- 
Kutta 4-th order method, similarly to work !4C?|. For 
computing the discrete Fourier transform, the FFTW li- 
brary was used J43j | . The length scale was normalized by 
a factor h = L/{2n) and the velocity scale by {ghy^^. 
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FIG. 3: (i) Free surface and bottom for t = 0. Only part of 
the entire periodic domain is shown. 
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FIG. 4: (i) t = 50: propagation stage. 



B. Numerical results for Eas. (l56ll - ll57l) 

The above remark about time-dependent w-period 
equally concerns the system (|^ - (|K7)l as (l5(^ - (|Hn|l . How- 
ever, for locaUzed disturbances the relative change of a{t) 
remains small. Therefore solitary waves are possible to 
simulate with a constant a and with the regularized T. 
It was done so the first numerical experiment [referred 
as (i)], where periodic boundary conditions were applied 
with the fixed u-period L = 200. The initial conditions 
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FIG. 5: (i) t = 69: breaking of the wave. 
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FIG. 6: (i) Crest of the wave for t = 69. 



were taken B{u^ 0) = and 



yl(M,0) 



z„(u, 0) 



2.5£;(u) + 0.25 1.4 



E(u)^\ C(w)' 



where 



E{u) = exp[8.0(cos(27rii/i)-0.1)], 
C[u) — cosh[(L/7r) sin{7r(u + i)/L}]. 

This initial configuration is symmetric (even). It results 
after some time in two oppositely propagating, nearly 
solitary waves. The waves are created in the region where 
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FIG. 7: (ii) The bottom profile and shape of the free surface 
at i = 0. The velocity field is everywhere zero. 



the depth is maximal. In the course of motion each wave 
approaches a shallow region where the surface profile y = 
ri{x,t) becomes steeper and finally multi-valued. 

Some of the results of this numerical experiment are 
presented in Figs l3l6l where the velocity distribution on 
the surface is shown, as well as shapes of free surface and 
of the bottom for several moments of time. In general, 
the computed wave profiles look quite realistic, though 
the present theory does not take into account viscous 
effects. Steeping of the wave profile is clearly seen. It 
should be noted, however, in these simulations the bot- 
tom shape is preserved only approximately due to dis- 
cretization errors, and the same concerns the total en- 
ergy. The computation was stopped well before the mo- 
ment of formation of a singularity on the crests of the 
waves, when the numerical scheme becomes invalid. In 
real world this moment corresponds to development of a 
three-dimensional instability resulting in vortices, splash 
and foam. 



C. Numerical results for Eas. (l59l) - II60B 

In the second numerical experiment [referred as (ii)], 
the shape of the bottom was fixed by analytical function 

2(0 = C + *^ exp(2^^(C + i)/L), 

with the dimensionless parameters L = 100, d — 0.16. 
The initial velocity field was taken zero: /3(u,0) — 0, 
while C(w, 0) had the form 

C(w,0) =7^ + 9(0.63, 0.0, u + i) + OM 9(0.9, 0.04, u + i), 

where 

9(r p w) ^ -^^ In ( ^ + ^ eM-^mjw - pL)/L) \ 

Qualitatively, these initial conditions are similar to those 
in the first experiment, however now two oppositely prop- 
agating waves are created over inclined region of the bot- 
tom, so there is no left-right symmetry in the evolution. 
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FIG. 8: (ii) Velocity distribution on free surface and shape of 
the surface at f = 6. 
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FIG. 9: (ii) The same as in Fig|Hl at t = 12. 
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FIG. 10: (ii) The same as in Fig|Hl at t = 17. 



The full system (|63|) - (|65|l was solved with a high accu- 
racy for TO in the limits —7000 < to < 7000 (the energy 
conservation was up to 12 digits for a long "smooth" ini- 
tial stage of the evolution, and it was up to 5 digits for 
the final stage just before the breaking). This more accu- 
rate numerical solution was compared to a less accurate 
solution obtained with fixed a — a(0). The difference 
was found very small. 

The corresponding numerical results are presented in 



Figs l7ll(?l Here we again observe steeping of the wave 
profile with the tendency towards finite time singular- 
ity formation on the crest. Such behavior indeed takes 
place in natural conditions when the flows are almost 
two-dimensional. 



V. SUMMARY 

In this paper we have derived approximate weakly- 
nonlinear, as well as exact nonlinear equations of motion 
for potential water waves over a strongly inhomogeneous 
bottom. The consideration was based on using the con- 
formal mappings. For linear waves over periodic seabed, 
the band structure of the spectrum has been calculated. 

Though the obtained exact equations can be written 
in formally the same form as those derived in [Sg for a 
straight horizontal bottom, but admissible solutions have 
different analytical properties if the bottom is inhomoge- 
neous. When the equations are written in this form, the 
bottom shape is preserved as an integral of motion. 

Numerical experiments have been carried out that con- 
firm advantage of the theory by giving quite realistic pic- 
tures for wave profiles before wave breaking. 

Of course, the above "inviscid theory" works only on 
large enough spatial scales and only until the singularity 
moment, as it was clear from the very beginning. Prac- 
tically, this theory is good for description sea and ocean 
waves before their breaking. 
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